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Abstract. A discrete Fourier analysis associated with translation lattices is 
developed recently by the authors. It permits two lattices, one determining the 
integral domain and the other determining the family of exponential functions. 
Possible choices of lattices are discussed in the case of lattices that tile R 2 and 
several new results on cubature and interpolation by trigonometric, as well as 
algebraic, polynomials are obtained. 



1. Introduction 

A framework of discrete Fourier analysis associated with translation tiling was 
developed recently in [7J, based on the principle that if SI is a bounded open set 
that tiles M. d with the lattice L — AZ d , then the family of exponentials {e 2ma ' x : 
a e L }, where L = A~ tr Z, d is the dual lattice of L, forms an orthonormal 
basis in L 2 (fl) ( 4 ). Our set up permits two lattices, one determining the integral 
domain and the other determining the exponentials that are orthogonal under the 
discrete inner product. The case that both lattices have the regular hexagon as 
fundamental domain was studied in [7J to illustrate the main set up, which leads 
to new cubature formula and Lagrange interpolation for trigonometric polynomials 
on hexagonal domains and equilateral triangles, as well as results for algebraic 
polynomials on the region bounded by Steiner's hypocycloid. This is extended to 
three dimension in [9] , giving results on cubature and interpolation on the rhombic 
dodecahedron and the tetrahedron, and further extended to M. d in [TU] for Ad type 
lattice. In [8], the two lattices are chosen differently with fundamental domains 
being a square and a rhombic (rotation of the square by 90°), respectively. The 
choice leads to, surprisingly, one family of minimal cubature for product Chebyshev 
weight on [— 1,1] , first discovered by working with common zeros of orthogonal 
polynomials of two variables. An extension to three dimension gives a family of 
cubature formulas on the cube that have the smallest number of nodes among all 
known formulas, which coincides, rather surprisingly, with the cubature discovered 
in [TT] by a totally different method. 

The two lattices in [8] were chosen for the purpose of obtaining algebraic cu- 
bature formulas on the square. Its success prompts us to ask what other choices 
are possible. In the present paper, we try to answer this question in the case of 
M 2 . Up to affine transforms, there are essentially two types of translation tiling 
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in R 2 with fundamental domain being either squares or regular hexagons. Their 
combinations in our discrete Fourier analysis, however, yield several distinct cases, 
including several cases not covered in our previous studies. One that is of particular 
interesting has one tiling sets as the regular hexagon and the other as the rotation 
of the regular hexagon by 90° (see Section 3.5), which leads to another set of cuba- 
ture and interpolation on the equilateral triangle, different from those obtained in 
[7] . In order to present the main idea without being overwhelmed by notations and 
numerous formulas, we shall work mostly with cubature formulas, a central part 
but by no means all of discrete Fourier analysis, unless other results are deemed 
novel enough to warrant inclusion. 

The paper is organized as follows. In the following section we recall the frame- 
work developed in [7] and use it to treat the classical product discrete Fourier anal- 
ysis on the plane, which illustrates well what can be expected in the non-classical 
settings. Section 3 is divided into a number of subsections, each deals with one 
distinct choice of two lattices. 

2. Discrete Fourier analysis with lattice 

In the first subsection, we give a succinct recount of the framework of discrete 
Fourier analysis with tiling in [7]. We refer to PQ for lattices, tiling and various 
related topics, and refer to [HQ3] for some applications of discrete Fourier analysis 
in several variables. In the second subsection, we illustrate the general theory by 
using it to recover the classical product discrete Fourier analysis on the square. 

2.1. Discrete Fourier analysis. A lattice L in R d is a discrete subgroup L = 
La ■= A7L d , where A, called a generator matrix, is nonsingular. A bounded set f2 
of R d , called the fundamental domain of L, is said to tile R d with the lattice L if 

Xn(x + a) = 1, for almost all x £ R d , 

where xn denotes the characteristic function of fl. We write this as + L = R . 
For a given lattice La, the dual lattice L\ is given by L\ = A~ tr Z d . According 
to a result of Fuglede [4] , a bounded open set £1 tiles R d with the lattice L if, and 
only, {e 27ria x ■ a 6 L- 1 } is an orthonormal basis with respect to the inner product 

(2-1) ( f) g) n = —}— f f(x)gjxjdx. 

mes(SZ) J Q 

For L = A% d , the measure of is equal to | det(A)|. Since L\ = A~ tx 7lf l , we can 
write a = A- tr k for «gL] and k E Z d , so that e 2 ™ ^ = e^^'A^x _ 

For our discrete Fourier analysis, the boundary of f2 matters. We shall fix an f2 
such that G and fl + AZ d = M. d holds pointwisely and without overlapping. 

Definition 2.1. Let £Ia and £Ib be the fundamental domains of AZ d and BTL d , 
respectively. Assume all entries of the matrix N := B tr A are integers. Define 

A N := {k e Z d : B^'k e Q A } and A^ := {k G Z d : A~ tr k £ n B }. 
Furthermore, define the finite dimensional subspace of exponential functions 
H N := span{e 2 ™ fc,rA ~ la: , k e A^} . 
The main result in the discrete Fourier analysis is the following theorem: 
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Theorem 2.2. Let A,B and N be as in Definition \2.1\ Define 

for f,g in C(TIa), the space of continuous functions on FIa- Then 

(2.2) (f,9)a A = (f,9) N , f,geH N . 

It follows readily that (|2.2[) gives a cubature formula exact for functions in Hn- 
Furthermore, it also implies a Lagrange interpolation by exponential functions. Let 
Tfqf denote the Fourier expansion of / in H.n with respect to the inner product 
(•, -)jv, which can be expressed as 

(2.3) I N f{x) = J2 f(B- tr kM B (x - B- tr k), f e C(Sl A ), 



feSA, 



where 



{ ' n * [ ' |dct(A0| ^ 

Theorem 2.3. Let A,B and N be as in Definition \2. 1\ Then Tn is the unique 
interpolation operator on N inTi^ ; that is, 

T N f(B- tr j) = /(B- tr j), Vj e A N . 

In particular, #Ajv = #A;ytr = | det(iV) | . The cubature formula and the La- 
grange interpolation are for functions that are periodic with respect to the lattice 
AZ d , which are functions satisfying 

f(x + Ak) = f{x) for all k e Z d . 
The function x i— > e 2T " fc A x is periodic with respect to the lattice AZ d . 

2.2. Classical discrete Fourier analysis. We deduce the classical result on the 
plane (cf. [2j[23]) from the general theory described above. As mentioned in the 
introduction, we shall limit our consideration to cubature formulas. The result 
hints at what is possible in the non-classical cases in the rest of the paper. 

For n G N, let A = I, the identity matrix, and B = 2nl. Then N = B tr A = 2nl 
has all integer entries. Let ft a — [— \, \) 2 ■, which tiles R 2 with Z 2 pointwisely and 
without overlapping. We shall write A„, Ajj, TL n in place of Ajv, A' N , TLn- Then 

A„ = A* = {k e Z 2 : k e [-n, n) 2 } and U n = spanje 2 ™' * : k E A*,}. 

It is clear that #A„ = (2n) 2 . The equation (|2.2|) in this setting becomes 

P -. n— 1 n—1 

(2.5) J i i 2 f(x)gjx)dx = — J2 E f(t>k)9(k>kl f>9^H n . 

[ 2 1 2 ] k\ =— n k-2 ——n 

To illustrate what can be done on cubature, we state the results in stages. 
Stage 1. It is easy to see that f|2.5p yields a cubature formula 

„ 1 n—1 

( 2 -6) / = 7T2 £ E f&>k)> 
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where 7i* := span{e 27rife :l; : k 6 [— n, n] 2 nZ 2 }. The set of nodes of cubature ()2.6|1 is 
not symmetric on [— |, ^] 2 , since it has points on only part of the boundary of the 
square. We would like to have cubature whose nodes is symmetric on the square. 

Stage 2. We construct cubature formulas for Ti^n-i that have symmetric nodes 
on the square. Such a cubature is invariant under sign changes in both variables. 
We can in fact obtain two such formulas from (|2.6p . The first one is obtained upon 
using the periodicity of the functions in the sums in the right hand side, 

n n 

l2 /(^ = ^2 E E <W(£>£). V/e^L-i, 

i 2-2l 2 fci=— nfc 2 =— n 

where Cfc >n = 1 if € (— n,n) 2 , cfc. n = 1/2 if either fci = ±n or &2 = ±n but not 
both, and Ck.n = 1/4 if fc = (±n, ±n). The second one is obtained by applying 
(|2.6p to the function /(• + -^) and using the periodicity of / in the integral, 

(2.8) f f(x)dx 

The fact that the set of nodes in either (|2.7|) or (12. 8j) is invariant under the group 
Z| (sign changes) allows us to derive cubature formulas for product cosine and 
produce sine functions. Let TC n :— spanjcos 2irkiX\ cos 2nk2X2 '■ < ki,k 2 < n} 
and TS n := spanjsin 2-nk\X\ sin 2^2X2 : 1 < fci, &2 < n}, which consist of functions 
in 7i* that are invariant or anti- invariant under Z 2 ,, respectively. 

Stage 3. Restricting (|2.8j) to TC^n-l, we obtain a trigonometric cubature, 

„ n—1 n—1 

(2-9) / = TT2 E E V / e ^n-x, 

whereas restricting ()2.7[) to TC2n-i gives another trigonometric cubature for TC2n-\ 
Furthermore, restricting (|2.7[) or (|2.8[) to TS^n-i leads to cubature for TSm-i- 

The Chebyshev polynomials of the first and the second kind are defined, respec- 
tively, by T„(£) = cos27m# and U n (£) — sin27r(n + 1)8/ sin27r#, where £ = cos27r# 
with < 9 < 1/2. These are orthogonal polynomials with respect to u?o(£) : = 
(1 — £ 2 ) -1 / 2 and wi(£) := (1 — S, 2 ) 1 ^ 2 , respectively, on [— 1, 1]. Consequently, under 
the change of variables 

(2.10) (x!,x 2 ) h-> (t/i,J/2) = (cos 2nxi, cos 2nx 2 ) £ [-1, l] 2 , 

the space TC n is mapped into the product space n„ x n n , where n„ denotes the 
space of algebraic polynomials of one variable, and TS n is mapped into n„_i xll n _i. 
Stage 4. Under the map x t—> y of (|2.10[) . the cubature (|2.9[) becomes 

, « -. n—1 n—1 

(2-11) ^2 / /(»)Wb(y)«^= i E E /("OB =^,006=^1), 

for / S n 2 „-i x n 2 n-i, where Wo(y) = ^0(2/1)^0(2/2), which is in fact the prod- 
uct Chebyshev-Gauss cubature of the first kind. Applying the same procedure 
on the cubature l|2.7[) . we obtain the product Chebyshev-Gauss-Lobatto cubature. 
Furthermore, if we apply this procedure on the cubature for TS^n—l that were 
mentioned in Stage 3, we obtain the product cubature for the product Chebyshev 
weight W\{y) = w\ (2/1)1^1(2/2) of the second kind. 



n—1 n—1 

= 4^2 E E ft*} 

k\——n ko=—n 
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3. Discrete Fourier analysis on planer domains 

We now apply the general theory in Section 2.1 on the non-classical choices 
of lattices. The guideline of our choices is the program that we outlined for the 
classical case in subsection 2.2. We list the cases according to the shapes of the 
fundamental domains of lattices. For example, the classical case in Section 2.2 is 
Square-Square. The main ones that we consider are the regular domains such as 
square, rhombus, and regular hexagon, which are depicted below. 



Figure 1. Square, rhombus, regular hexagon 

3.1. Square- Rhombus. In this case we choose A = I with £Ia = [—5, \) 2 being 
the square and choose B = nR, where i?Z 2 has rhombic as its fundamental domain, 

with Qb = {x £ R 2 : — n < x\ + x 2 < n, —n < x 2 — x\ < n}. 

This case was studied in [8]. We shall be brief. Here = Ajy =: A n , where 

A„ = {j S Z 2 : -n < 32 ± h < n} and A* n := {j £ Z 2 : -n < j 2 ± ji < n}. 

The set A„ is not symmetric on [— n, n] 2 but A* is. The cardinality of A n is 
|A„| = 2n 2 . We follow the program in Section 2.2: In Stage 1 we deduce a cubature 
from Theorem 12.21 which has nodes indexed by A„, then in Stage 2 we derive a 
cubature by periodicity that has nodes indexed by A*. By considering functions 
that are even in both variables, we deduce in Stage 3 a cubature for trigonometric 
polynomials, which we state as follows. Changing variables from j to k = 2nB~ tr j, 
or fci = ji + ]2 and k 2 — J2 — ji , it follows easily that j £ A* is equivalent to 

k £ X* = {2k : -f < ki,k 2 < f}U{2fc + l : -2±1 < ki,k 2 < ^} . 

Let X°, X!^ and X% denote the set of points in X* that lie in the interior, the edges 
excluding corners, and the corners of [— n,n} 2 , respectively. 

Throughout the rest of the paper, we will adopt the convention that X°, X e and 
X v are subsets of X defined as above, whenever the domain to which the interior, 
edges and corners relate to is clear. 

Theorem 3.1. For n > 2, the cubature formula 

r 1 fl, k£X° 

(3.1) / f(x)dx = _^ C W/(A) with c^= I, k£Xl 

is exact for f £ T 2n _ x , where T* n := span { e 27rlk x ■ k £ A* s } . 



R := 



1 1 

-1 1 
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The index set A* is most suitable for dealing with algebraic polynomials. In 
fact, under the change of variables x i— > y in (|2.10j) . the space TC* m becomes the 
space II^j of algebraic polynomial of total degree m and the cubature (|3.ip becomes 
a cubature for the product Chebyshev weight Wo(y) that is exact for n| n _ 1 . Let 

H„ := {(2fci, 2fc 2 ) : < k u fei < f } U {(2fci + 1, 2fc 2 + 1) : < k u k 2 < 2=1}. 

Then, in Stage 4, (|3.1[) becomes the following: 

Theorem 3.2. Le£ = (cos cos ^ ZL ). ITiera i/ie cubature below is exact for 

(3-2) - / /(y)Wo(iOdy = J] A fe n) /(^-), := \ 2, fee S«, 
^ • /[ " 1,1]2 " fc es„ [i, fcesj. 

The cardinality of 3„ is |S n | = "^" 2 +1 - > + |_§j + 1, which is just one more than 
the theoretic lower bound for all such cubature ([3] [12]). The formula (|3 . 2[) first 
appeared in |20j . where it is constructed by considering the common zeros of or- 
thogonal polynomials of two variables; see also [13]. We can also derive similarly 
cubature for the product Chebyshev weight W\(y) of the second kind. 

The Lagrange interpolation polynomials based on the points in S„ were defined 
and studied in [3T]. The result there has also been recovered in [5], with (|2.3p as 
the starting point, by utilizing the discrete Fourier analysis. 

3.2. Rhombic-Square. In this case we choose A = R with fundamental domain 
&r = {x G K 2 : — 1 < X2 ± x\ < 1}, the rhombic, and B = nl. Again write A„ ... 
in place of Ajy It is then easy to verify that A„ = with 

An = {j G Z 2 : -n < k 2 ± fci < n} and A; := {j G Z 2 : -n <k 2 ±k 1 < n}. 
Furthermore, the space of exponential functions Ti. n is given by 
H n = {e^i+Mzi+fe-fci)^) . k e An j. 

and H* n is likewise defined in terms of A* . Changing variables shows that 
H* := { c ^U^+^2) . _ n < jl h < nj h = n (mod 2 )}. 

Following the program in Section 2.1, it is easy to see that the cubature in Stage 2 
that has symmetric set of nodes, indexed by A*, takes the form 

i r i f 1 ' fceA "' 

(3.3) - / f(x)dx = — 2 Y. c ^f(±), f e W^_ x , <%,„ = i fc G A£, 

The subspace of functions in 7i* that are even in both variables becomes 

(3.4) T n := span{cos7rjiXi cosirj 2 x 2 : < ji, j 2 < h = h ( m od 2)}. 

For functions in 72„_i, we only need to consider the triangle Tr := {x : < 
xi,X2, xi + x 2 < 1}. Thus, in Stage 3, cubature (|3.3p becomes 

(3-5) 2/ /(x)d!c = ^ Yl **.«/(£)> feT 2n _ r , 

where Tr is the triangular domain S„ = {(fci,fe) : < k\,k 2 ,k\ + k 2 < n}, and 
Afc, n = 4 if fe G H°, Afc,„ = 2 if k G E£, A(o,o),n = 1) an d A( nj o),n = A( ,„),„ = 1/2. 
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Under the mapping x > y = (coS7ra;i, cos7nE2), the boundary x\ + X2 = 1 
of the triangle Tr is mapped onto yi + y 2 = 0, so that Tr is mapped onto the 
triangle Tg = {(y^y 2 ) : — 1 < 2/1)2/2 < 1, 3/1 +2/2 > 0}, which is half of the square 
[— 1, l] 2 . The cubature (|3.5|) in Stage 4 becomes a cubature with respect to the 
product Chebyshev weight Wo over Tg that is exact for the subspace of polynomials 
ni n -i = {T fcl (a: 1 )T fc!1 (a;a) : < fc i> fc 2 < 2n - l,ki = k 2 (mod 2)}, the image of 
T-in-i under the same mapping. Since II* does not contain polynomials of total 
degree, we shall not write this cubature explicitly out. It is easy to see, however, that 
this cubature is in fact half of the product Chebyshev-Gaussian-Lobatto cubature, 
in the sense that its domain is half and it is exact for half of the polynomials of the 
latter cubature. 



3.3. Rhombic-Rhombic. Here we choose A = R and B = nR~ tr — ^A, so that 
N = B tr A = nl have integer entries. Then f2^ = Or as in the previous case. Again 
denote A^v, ... by A„, ... . It is easy to see that A n = AJ, with 

A„ - {] G I? : -f < -ji,h < f } and A; := {j G Z 2 : -f < n ,j 2 < §}. 
Moreover, the space of exponential functions TL n is given by, as in Section 3.2, 

H n — { e Tr i((kl+k 2 )x 1 + (k 2 -k 1 )x 2 ) . y, g ^\ 

and is likewise defined with A„ replaced by A*. In this case, the cubature 
derived from Theorem l2.2i in Stage 1, takes the form 

(3-6) ~ J f{x)dx = 1 J] /(^, V/ G H* 2n _ v 



keA„ 



The set of nodes of this cubature is on Or, and it contains no points on the boundary 
of fifl; when n is an odd integer, whereas it contains points on half of the boundary 
of flu when n is an even integer. In the latter case, we can again derive a cubature, 
exact for W|n-i) tna * nas n °tes indexed by A* as in Stage 2. Let us consider, 
however, only the case of n being an odd integer below. As can be seen upon 
changing variables ji = k\ + k 2 and j 2 — k 2 — k\ , the subspace of functions in W* 
that are even in both variables is exactly T n in (|3 .4(1 . Thus, just like in the case 
of Rhombic-Square, restricting ()3.6() to functions in Ti. 2n _ 1 that are even in both 
variables, we deduce a cubature of Stage 3 on the triangle Tr, 



!•'--) 2 [_ f(x)dx = ^Y1 ^nf(^, V/ G r 2n _ ll 

fcGH„ 



where H„ = {(ki,k 2 ) : < hi, k 2 , fei + k 2 < f }; A fei „ = 4 if k G 3°, A fcj „ = 2 if 
either k\ = or fc 2 = or fcx + k 2 = 5 but not both (i.e.,fc G 3^), and A(o,o),n = lj 

\o,n),n = 1 • 

Finally, under the mapping x ^ y = (cos irx\, cos 772:2), the cubature (|3.7p be- 
comes a cubature with respect to the product Chebyshev weight Wo over the trian- 
gle domain Ts for the polynomial subspace Il2 n _ 1 defined in the previous subsec- 
tion. This cubature is exactly half of the algebraic cubature in the Square-Rhombic 
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3.4. Hexagon-Hexagon. In this case we choose A = H and B — %H, where 

H= {^1 2) With n H = {xeM. 2 : -l<x 2 ,^f±±f <l}. 

This case was studied in [7J . We shall be brief, but recall necessary definitions that 
are needed in the following subsection. As shown in [7j [19] , it is more convenient 
to use homogeneous coordinates (ti,t 2 ,t 3 ) defined by 





(4 


ii) ■ 1 







\-4 




which satisfies t\ + t 2 + t 3 = 0. We adopt the convention of using bold letters, such 
as t to denote points in homogeneous coordinates. We define by 

R 3 H := {t = (ti,t 2 ,t 3 ) e W 3 : ti + 1 2 + <3 = 0} and M:=Z 3 r\M. 3 H 

the spaces of points and integers in homogeneous coordinates, respectively. In such 
coordinates, the hexagon becomes 

fl = {t £ R 3 H : -1 < t u t 2 , -t 3 < 1} , 

which is the intersection of the plane t\ + t 2 + £3 =0 with the cube [— 1, l] 3 . The 
index sets A„ and Ajj satisfy A n — A' n = H„, where 

H„ := {j e H : -n < j u j 2 , -33 < n} and H* := {j e H* : -n < j uj2 , -j 3 < n}. 

Furthermore, since, for k = (fci, k 2 ), k tr H~ 1 x — -|k-t with k = (k\, k 2 , — k\ — k 2 ) tr E 
H, the exponential functions and the space 7ijy become 

0k(t) := e^ k,rt and H n := {0 k : k e H„}. 

In homogeneous coordinates, x = y (mod H) becomes t = s (mod 3), which is 
defined by t\ — s\ = t 2 — s 2 = t% — S3 (mod 3), so that / periodic in H, i.e. 
f(x + H) = f(x), becomes /(t) = f(t + j) whenever j = (mod 3). 

In this case, the cubature derived from Theorem 12.21 in Stage 1 has nodes over 
{^:kG H n }, from which we derive another cubature, the set of nodes of which is 
symmetric and indexed by H* , as in Stage 2: 

Theorem 3.3. The following cubature is exact for f € 7~L 2n _i, 

The group of isometries of the hexagon lattice is generated by the reflections 
in the edges of the equilateral triangles inside the regular hexagon, which is the 
reflection group A 2 . By considering the invariant and anti-invariant functions under 
A 2 in the space TL ni we end up with functions that are analogues of cosine and sine 
functions on an equilateral triangle, and the cubature (|3.9|) becomes a cubature on 
the triangle for such functions. To be more precise, we choose the triangle as 

(3.10) A :={(i 1( t 2 , t 3 ) : t x + t 2 + t 3 = 0, < t x , t 2 , -t 3 < !}■ 

The region A and its relative position in the hexagon are depicted in Figure [21 
where the points are labeled in homogeneous coordinates. The generalized cosine, 
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Figure 2. The fundamental triangle of the regular hexagon. 



TCk, and the generalized sine, TSk, are defined in terms of 

(3.11) V+f(t) := ]T f(ta) and V~f{t) := £ sign^ta) 

ueA 2 &eA 2 

as TCk(t) := V + 4'w{'t) and TSk(t) := V~^>k(t), respectively; more explicitly, 
1 



(3.12) 



TCk(t) 



3 



e s 



(k 2 -k 3 )(t 2 -t 3 ) 



cos k\-nt\ + c 3 



+ eT( fe 2-fe 3 )(*i-t2) C0S k lW t 3 



(fe2-fe 3 )(t3-tl) 



ke A, 



cos k\Txt 2 



(3.13) TSk(t) =- 



e ^(*2-fc 3 )(t a -t 3 ) sin fc l7ril + e ^(fa-fcs)(t3-ti) S i n fc l7 r< 2 



-e 3 



(fc 2 -fc 3 )(ti-*2) 



sin fci7r<3 



ke A°, 



where A := {k e H : fci > 0, /c 2 > 0, fc 3 < 0} and A° is the interior of A. These 
functions are orthogonal with respect to the integral over A, and they are elements 
of TL n that invariant and anti- invariant under A 2 , respectively. The cubature Q3.9P 
when restrict to invariant functions becomes, as in Stage 3, the following: 

Theorem 3.4. Let TC n := span{TCk : k € A n }. The cubature below is exact for 
all f e TC 2n _i, 

, 1 n h f6, jGA°, 

(3.14) 2 / /(i 1 ,f 2 ) Wa = - ? ^^Af ) /(if), Af := 3, 

*=°*= i, jeA;. 



The nodes of the cubature (|3 . 14[1 are equally spaced points in A (Figure 2). 
The generalized cosine and sine functions can be mapped into algebraic polyno- 
mials of two variables under the following mapping, 



(3.15) 



cos f (t 2 - h) cos | (t 3 - t 2 ) cos | (ti - i 3 ) - ±. 



a; — 3 

y = fsinf (t 2 - h) sinf (i 3 - £ 2 )sinf(ii - t 3 ), 



which are the real and imaginary part of TCo,i,-i (t) , the first non trivial generalized 
cosine function. Under this mapping, we call the polynomials 



T^(x,y) := TC fc , m _ fc ,_ m (t) and UJ?(x,y) := 



TS 



fc+l,m— — m— 2 

TS M ,_ 2 (t) 



(t) 
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where < k < m, generalized Chebyshev polynomials of the first and the second 
kind, respectively. They are algebraic polynomials of total degree n and are or- 
thogonal polynomials with respect to the weight function w_i(x 1 y) and 
respectively, where w a (x,y) is defined by 

w a (x, y) = ^tt 4q [~3(a; 2 + y 2 + l) 2 + 8(x 3 - 3xy 2 ) + 4]" , 

and the integral domain is the region A* bounded by the Steiner's hypocycloid, 
depicted in Figure 3, which is the region on which w a (x,y) is positive. These 




Figure 3. The region A* bounded by Steiner's hypocycloid. 

polynomials were first studied in [6|. As in Stage 4, the cubature (|3.14j) under the 
change of variable (|3.15[) becomes a cubature for w_i(x,y)dxdy on A* that has 
dim II 2 nodes and is exact for algebraic polynomials of degree 2n — 1. Furthermore, 
we can derive a cubature from (|3.9p for anti- invariant functions in Stage 3, which 
becomes under (|3.15p a cubature for wi(x,y)dxdy that has dimlT 2 ^ nodes and 
is exact for algebraic polynomials of degree In — 1. The latter one provides an 
example of a family of Gaussian cubature formulas, a rarity of only the second 
example known so far (the first one appeared in [17]); see [7] for details. We refer 
to [31 [TBI [18] for the topic of Gaussian cubature. 

We now address one question that was not addressed in [7] . Taking the cue form 
the cubature 12.81 in the Square-Square case, we can apply the cubature derived in 
Stage 1 on the functions fit + a) and then use the hexagonal periodicity of the 
integral to derive the following cubature in Stage 2, 



(3.16) 




and hope to choose a so that the set of nodes in (|3 . 1 6[) is symmetric. The question 
is if it is possible to find a a so that the set of nodes has full symmetry of A2 ■ 

It is easy to see that if a = (01, et2, — a% — 0,2) £ R# satisfies \a\ |, |aa| < 1/n, then 
the set of nodes of (|3.16|) will be inside the hexagon f2, although not symmetric for 
most of the choices. The two cases that offer the most symmetry are 

ai := (^T'^r' _ Jr) and a 2 : = ( _ ^r>-^r> 

where, when a.2 is used, we need to use the periodicity of / (or congruent relation 
with respect to H) to make sure that all points in (|3.16p are in fL Neither of these 
two choices, however, offer complete symmetry under the group Ai- In Figure 4, 
we depict the set of points resulted from these two choices. Each set of the points is 
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(-n,0,») 




FIGURE 4. Nodes of (pUB]) with a x (left) and a 2 (right) 

invariant under a subgroup of A 2 of three rotations, but neither is invariant under 
the group A 2 . As a result, we cannot restrict the cubature ()3.16|) with either ai or 
a 2 to the generalized cosine or sine functions in hopes of obtaining new cubature 
on the triangle in Stage 3, in contrast to Square-Square case. 

The interpolation on the hexagon and on the triangle were studied in [7j . In par- 
ticular, we have Lagrange interpolation based on equally space points on the triangle 
A, which enjoys a closed formula in trigonometric functions and has Lebesgue con- 
stant in the order of (logn) 2 . One can also consider approximation on the hexagon 
and the triangle ([H]) for functions that are periodic in H . 

3.5. Hexagon-Hexagon Transpose. Here we choose A — H, the matrix for the 
hexagon lattice, and choose B = nH~ tr with n E Z, so that N = B tr A = nl has 
all integer entries. The fundamental domain of the lattice Lb is given by 

nB = {xeR*:-^<x 1 ,fx 2 ±lx 1 <^}. 

3.5.1. Cubature. It is again convenient to use homogeneous coordinates as defined in 
the previous subsection. The fis is the regular hexagon in Fig. 1 rotated by 90°, as 
depicted in Figure 2, in which the right hand figure is labeled in homogeneous coor- 
dinates. Here the index set Ajy — {k 6 Z 2 : — n < Iki + k\, k-x — k\, 2k\ + k 2 < n\, 




Figure 5. The hexagon domain and Sis 



which becomes in homogeneous coordinates defined by 

Kl := {j G H : -n < 32 ~ ji,ji - ja,j2 - h < n}. 

We also have Ajv = {k S Z 2 : — n < 2k 2 — k\, k\ + k 2 , 2fci — k 2 < n}. Recall that 
t = (mod 3) means, by definition, t\ = t 2 = t% (mod 3). It is not hard to see 
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that the set A at becomes, in homogeneous coordinates, K n defined by 
Kn:={j£B:-n< j u j 2 , -j 3 < n, j = (mod 3)}. 

We also denote by Kjj* and K* the sets defined with < in place of < in and K„, 
respectively. The set K* can be obtained form a rotation of K£* , as shown in the 
following proposition, which can be easily verified. 

Proposition 3.5. For k = (fci, k 2 , G H, define k := (ka — k 2 , ki — k%, k 2 — k\). 
Then § G i/ke K£ and k G K* i/kG Kj* . 

The finite dimensional space TL n of exponential functions becomes 

K n := span j<£j(t) = e^"* : j G } and /C* := span { 0j : j G IK**} . 

By induction, it follows that dim/C„ = |K n | = n 2 and dim/C* = |K* | = n 2 + n + 1 
if n = 0, 2 (mod 3) and |K* | = n 2 + n — 1 if n = 1 (mod 3). The two sets K% n and 




Figure 6. The set K% n (left) and the set K f 3 * n (right). 
K3* take different shapes, which we depict in Figure [6l Define 

(3-i7) (/,</>« : =^E /(*)7(J) and (/.ff>; = 4 £ tf°/(£)7(J)> 

jeK„ jeK* 
where Cj are as defined in (|3.9p with H„ replaced by K„. 

Theorem 3.6. for » > 0, |^ J n f(x)g(x)dx = (f,g) n = (f,g)* n for f,g G /C„. 

In particular, (^>j,0k)n — 1 */j = ^ (mod 3n) and (^jj^k)^ = otherwise, for 
k,j G K„. Moreover, we have the cubature 



(3.18) 




!n-l- 



The part of the theorem on (/,<?) n is exactly Theorem 12.21 while the part on 
(/, <?)* and the cubature can be proved by periodicity, just like the proof of Theorem 
13.31 in T| , upon using the Proposition 13.51 The cubature (|3 . 18[) is already one in 
Stage 2; we can also derive a cubature with nodes indexed by K„ as in Stage 1. 

Next we consider the invariant and anti-invariant functions under A 2 , which 
are the generalized cosines TCk and the generalized sines TSk considered in the 
previous subsection. By restricting to such functions, we again obtain cubature on 
the triangle A. The index set of the nodes of the cubature, denoted by T, is 

T„ := {j G H : < j u j 2 , -33 < n, j = (mod 3)} 
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derived by symmetry from K* , whereas the index set of the invariant functions 
being integrated exactly by the cubature, denoted by TT, is derived from Kj,*, 

Tj, = {j e H : < j t , j 2 , -n < n ,h - to < n,ji - j 3 < n}, 
which is inside a quadrilateral; Figure [7] shows its relative position in f^. We define 




V 3 ' 3 ' 3 ' 

Figure 7. The fundamental quadrilateral of under A2 



the following subspaces of trigonometric functions, 

TC n = span{TC k : k e Tjj and TS n = span{TS k : k 6 T^ }. 

The set T„ takes a symmetric form when n is a multiple of 3. In Figure ([8]) we 
depict the index sets T3„ and Tg n . 

(0,3n,-3n) Ifi.^f-'-f) 



(0,0,0) (3n,0,-3n) ( - ' > (**■*) 

Figure 8. the index sets T 3 „ (left) and T\ n (right). 



Theorem 3.7. The following cubature is exact for all f G TC2n—i> 

1 r 1 f6, J6T°, 

(3.19) — / /(t)A = - 2 2 Ai n) /(^), Af) := 3, j e T^, 

The formula (|3.19p is derived from p,18p by using the invariance of the functions 
in TC 2 „-i and the fact Q = ( U ae A 2 {ta : t e A°}) U ( U a eA 2 {ta : t e OA}). As 
the proof is similar to that of (|3.14[) in |7j, we shall omit the details. 

Similarly, we can also derive a cubature for TS 2 n -i based on points in T° . These 
are cubature in Stage 3. We note that the set of nodes in (|3 . 19|) is different from 
that in (|3.14p . see Figure 2 and Figure 8, even though both are on the triangle. 

As in the case of Hexagon-Hexagon, we can continue to Stage 4, where the 
cubature (|3. 19|) is mapped by the change of variables (|3. 15(1 to an algebraic cubature 
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for W-i/2{x,y)dxdy on A*, the region bounded by Steiners hypocycloid, which is 
a cubature exact for all polynomials in n| n _i but with many more nodes than the 
one derived in the Hexagon-Hexagon case. 

The set of cubature points in (|3.19[) and its image in the region bounded by 
Steiner's hypocycloid have also been studied in [15] . 



3.5.2. Interpolation. Applying Theorem 12.31 to the current set up, we obtain an 
interpolation operator X n f that interpolates / on K„ in the hexagon. We would like 
to consider interpolation on the triangle based on points in T„. For this purpose, 
we first construct a near interpolation operator on the symmetric set of points K* . 

Theorem 3.8. Let Sj := {k E K* n : k = j (mod 3n)} . For f e C(Q), define 

Kf(t) ■= E /(i)*»(t - i). *«(*) 4 n) ^ (t) - 

Thenl*J e K* and 2£/(i) = /(I) if j e K°, and 2*/(i) = £ keSj /(|) 
j € 9K* , £/ie boundary of K*. Furthermore, <&* (t) is a real function and it is given 
by the following formula when n = (mod 3), 



(3.20) $„(t) =4 



(2 cos 7TSt + cos nti) sin Trtj cos 27r ^' 
— ' sin nti sin 7rf2 sin irt^ 

1 



1 3 

IE 



i=l " 

cos + cos + C0S 



3 

w/iere si = ^^,52 = h i^,s 3 = h ^ ■ 

Proof. By Proposition 13.51 k S K* implies | G Kj*, and j e Kj,* implies j G K*. 
By homogeneity, k • t = (£3 — &2)si + (fci — fcs)s2 + (&2 — &i)s3 = k • s. As a result, 

**(*) = \ E 4 Tl) ^( s ) = -4 E c i"Vj(s)- 

?i * — ' J J n *■ — ' J 

Consequently, by the definition in (|3.17p . 

*«(¥) = ^E 4"Vi(¥) = ^E <h (B Vj(£) = te,^. 

IGK* IGK* 



Thus, by Theorem 13.61 it follows that 



(3.21) d>„(^) = <^,<%> 



0, otherwise, 



which proves the stated result of I*/. 

To derive the compact formula for we essentially need a formula for the 
Dirichlet kernel, denoted by 9„(t), of the Fourier series over H^*, 

e n (t) = E <PM 

Indeed, by the definition of c|j , it follows that 

!/^^,^ ^ fTCn,a,_ap(t), n = (mod 3), 



(3.22) n 2 $„(t) = -(0„(t) + n _i(t)) - 



3 1 3 ' 3 



2 10, otherwise. 
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Using the identity k • t = k • s and Proposition 13. 5[ we derive that 

e»(t) = £ Mt) = £ £ (s) = ]T 

We now partition K* into three parts according to the congruence relation, 

14°) := {j € K; : ji = j 2 = J3 = (mod 3)}, 
K« := {j S K* : ji = j 2 = j 3 = l (mod 3) } , 
< 2) := {j S K* : j% = j a s j 3 = 2 (mod 3) } . 

Using the fact that j e Ki "* | S H*„j, where H* is the index defined in the 

previous subsection, and 0j(s) = (3s) if j € K„ , we obtain from the Dirichlet 
kernel over H* in (3.10) of [7J, 

1 S i n7r |_»±3j^ 3 sinTrLfJ^- 



E*w=n^F-n 



sm7ri,- - LJ - sin7rt, 

icK (») 3 = 1 J 3=1 J 

(3.23) jeK " 

_ 3 sin27rLf Jtj - sm2 7 rL 2 ^Jtj 
r - ^ 4shi7rfi sin7rf 2 shi7rf3 ' 

3=1 

where we have used the identity sin 2^ + sin2i 2 + sin2t 3 = — 4 sin fx sinf 2 sini 3 
in the last equal sign [7J (3.15)]. Next we note that can be divided into the 
following three (non-overlapping) subsets = I&l 1 ' 3 '' U K^ 1 ' 2 ^ U Ki 1 ' 1 ', where 

K n 4) = {(~h - js.ja.Js) : h = 33 = 1 (mod 3), 1 < j 2 , -j 3 < n} , 
4 1,2) = {(jl.-ji - 33,3s) ■ 33 =ji = l (mod 3), 1 < j 3 , -ji < n} , 
Kl 1,3) = {O'l.ia, -ji - j 2 ) : ji = j 2 = 1 (mod 3), 1 < j u -j 2 < n} . 



Using the last set K„' 3 ', we define 



I(ii,t a ):= £ 0j(b)= £ £ 



e T ZL (32tl-3l*2) 



J 71 3|ji-1 3| j 2 — 1 



(3.24) _ e-^q - e- 2< * W«') e -^(l - e - 2 ™ W'*) 

1 g— 2i7r*2 ^ _ g— 2i7rti 

= (1 - C-^Oq - C- 2 "*2)(l - g-2^t 3 ) ' 

where the second equal sign follows from j-s = ji(si — S3)+j 2 (s 2 — S3) = j 2 ti~ ji^2- 
Moreover, we have j • s = fah - j' 1 t 2 = js*2 - j2*3 = ji*3 - j3*i. which yields 

T(t 2 ,t 3 )= ^ and I (*3,ii)= £ ^( s )' 

As a result, we conclude 

£ <f> 3 {s) =l(tx,t 2 ) +l(t 2 ,h) +l(t 3 ,ti). 
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Furthermore, we note that Kn* 1 = |— j : j G and, consequently, 

£ 0j(s)= X! <^ s ) + E <^( s ) = E <^( s ) + E ^j(- s ) 

jeK^UK™ jeK^ 1 ' jeK<?> jeKSj 1 ' jG^ 1 ' 



£ 0j(s)+ E 0j(s) = 2K{2'(ti,t2)+2-(t2,*3)+2'(t3,*i)}. 



Now assume that n is a multiple of 3. By using (|3.24|) and the fact that t is 
homogeneous, we obtain 

^ 0j(s) = 2R{I(t 1 , ta)) + 2R(J(t 2 ,t3)) + 2»(T(t 3) ti)} 



= 23? 



2i7rs 2 _ g— 2iirai\/j _|_ g2i7r-g-t 3 _ g- 





(1 


- e~ 2i ^){l ■ 


_ e -2»7rt 2 )(l 


_ g— 2i7r*3\ 


(e- 


-2Z7TS3 


e - 2i7 " 2 )(l + 


e 2l7r f 4l - c~ 


-2i7T^-t 2 _ g-2i7T^-t 3 ^ 




(1 


- e~ 2i7rt i)(l 


- e" 2i7rt2 )(l 


_ g— 2iirt 3 ^ 


(e- 


-2iiTSi 


e -2i7T S3 )(l + 


g2i7r-g-t 2 g— 


2iw%t 3 _ g— 2i7T-|ti^ ' 



(1 - c- 247rt i)(l - e^ 247r *2)(l - e" 2 ^) 

Combining the numerators and collecting the terms in 1, e 2l7T t' and e _2l7r t', we 
obtain that the combined numerator is equal to 

_ g-2i7TS2^ jg-2i7rsi _ g-2i7rs 3 \ 

/^g— 2J7TS3 _ g-2i7TS 2 ^ _|_ ^g-2l7TSi _ g— 2i7TS3^g— 2l7T Jt 3 

/^g— 2i7rsi _ g— 2j7TS3^ _j_ ^g-2i7rs 2 _ g — 2j7rsi 2i7r§ti 
^ e -2i7T« 3 _ g-2i7r Sl ^ + ^g-2i7rs 3 _ g-2iirs 2 ^\g-2i7rf t a 
3 -2i7r-Jt 3 -j _|_ /g-2t7rs 3 _ g-2i7rs 2 \ fg2i7r§ {j _|_ g2i7r^ti^ 



^g-2i7TS 2 


e" 


-2i7rsi^ _|_ /g— 2i"7rs 3 


■ ( C - 2 " S2 




g-2i7TSi ^g2i7T^-t 3 _ 


+ ( e - 2i7rs3 




g-2i7rs 2 \g2i7r-^ti _ 


+ (c- 2 ""! 




g-2i7TS 3 -jg2i7T^t 2 _ 


= (e _2i7rS2 — 


e" 


-2i7rsi\/g2i7rft 3 + 






g — 2l7TS3Wg2i7T-jt2 


= 2 cos 2w ™ t3 


gI7TS3 /gWrt 3 _ g— i-rrt 3 



e -2i7Tf t 2 ~j 



2 cos 2™*i e i7r;il (e i7rtl - e"™* 1 ) 



2 cos 2zl^*2 e " S2 (e i7r * 2 - e"' i7r * 2 ) 



3 

= 4ie l7rS3 cos s in7r<3 + 4ie iTrsi cos siliTrti + 4ie l7rS2 cos 2 -^ shi7rf 2 , 

where we use the facts that £3 = si — S2 , t\ = S2 — S3, £2 = S3 — si and s\ + S2 + S3 = 
for the second equal sign. Using £1 + £2 + £3 = 0, the denominator becomes 

(l_e- 2 ^*i)(l _ e - 2i,r * 2 )(l - e ~ 2M3 ) = -8isin7rtism7rt 2 sin7r£ 3 . 

Consequently, we derive that 

E <M s ) 

jeK^'uK*, 2 ' 

ins 3 2™*, in7rf + iv3l Zjmti • t + e irs 2 2™^ j f 

— 23? 3 - 2 - 2 1 

— 2 sin 7r£i sin 7r£2 sin -nt^ 

cos — g-^- cos irsj sm 7rr j 
^-^ sin 7rfi sin 7r£ 2 sin 7r£3 
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Combining the above equation with ()3.22j) and (|3.23|) . we obtain 



2 $„(t) = -(e^ + e^it)) -Tc f)f> _ ¥ (t) 

E <^ s ) + E ^(s))-Tc SiSi _ ¥ (t) 



27r(n— 3)ij 



sin 



27r(K+3)tj 



3 27rnti . , 

cos — g-^- cos 7TSj sin irtj 



8 sin Trti sin 7rf 2 sin 7rf 3 



E 



sin 7rtx sin 7rf 2 sin 7rf 3 



TO 



.(t) 



3 2lT7lt-i . , . 

COS — 3-^ Sill 7TEj COS 7ITj 



E 



3 27tnt~ . , 

cos — 3-^- cos irsj sin 7rtj 



2 sin 7rfi sin 7rf2 sin 7rf 3 



EWE 



i'=l 



sin 7rf 1 sin 7rf2 sin 7rf 3 



^ cos 2 ™ tj sin ntj (cos 7rtj + 2 cos nsj ) 1 . 27rn{ 
^— ^ 2 sin 7r^i sin 71^2 sin ^3 3 



TC, 



27rn^2 



>(t) 



27TKt 3 ' 



2 sin 7r£i sin 7rf2 sin 7rf 3 
This completes the proof. 



COS + cos — 3 



□ 



We now proceed to interpolation on the triangle A. The idea is to use the 
periodicity and apply the operator V ± in p. lip on the interpolation I n f, as in 
Theorem 4.7 in [7]. First we apply V~ on T n f , which gives the following: 



Theorem 3.9. For n>0 and f € C(A) define 



£ n /(t) := £ /(^j,„(t), 



$,»(*) = £ E AL n) TS k (t)TS k (i), 



kexj. 



where 



6, 


ki,k 2 ,n + 


fc 3 - fci 


n + fc 3 - fea > 0, 


1, 


k = 0, 






< 2, 


\e — ( n n 

V 3 1 3 1 


2n\ 
3 /' 




3 

2 ' 


k = (|,0, 


-f) or 


CO 2 -S) 


[3, 


otherwise 







T/ien £ ra is i/ie unique function in TS n that satisfies C n f(j-) = f(^), j € T° . 
Proof. By the definition of P and TSk, 

^(t-^) = i E 4 3n) ^(i)n-^(t) = A 53 „<*o- • 

keKj 

M 2 | 



„_ _ ' k *k(*)TS k (t) 



(3n); 



4 E E \ 

keX t* a£A 2 



w(i)TSka(t) 



I-A2 



£ c ( 3 ">M 2 |TS k (t)-L £ signCa)^) 



rr ' k 



1 E 4 3n) l k ^l TS ^ t ) TS ^^) = ^ E A k " ) TS k (t)TS k (i). 



kexj, 



kexj;* 
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Now, forj,leT°, 



= = ± 

crG^ 2 keK^* 3 3 aeA 2 iGK* 3 3 

= \ E = J E w - J^i. 

o-e^ 2 3 3 o-eAj 

where (Jj^ equals 1 if £ = J£ (mod 3), and is otherwise. This completes the 
proof. □ 

In fact, f? n (t) = 6'P _ $ n (t — i), where "P" acts on the variable t, from which 
the proof reduces to verify formula of £? given in the theorem, using the periodicity 
and the symmetry. Applying now V + to T n f, we obtain similarly the trigonometric 
interpolation on T„ in A. 



Theorem 3.10. For n > and f e C(A) defir 



A («) 

C/(t) := E /(iXUt), ^j.n(t) ^ E A k n) TC k (t)TC k (i), 
je T « kexj, 

where Aj are defined in (|3.19|) . T/ien £* is £/ie unique function in TC n that 
satisfies £*/(!) = /(i), j € T n . 

For n being a multiple of 3, we can deduce a compact formula for £? n (t) and 
^j, n (t) from that of ()3.20p . The interpolation points of £*/ are depicted in Figure 
8. From the explicit formula of <&* in (|3 . 20(1 . it is not difficult to prove, following 
proof of Theorem 3.6 in J7J, that the uniform operator norm (Lebesgue constant) of 
T*f in Theorem 13.81 satisfies ||Z*||oo < c(logn) 2 for n = (mod 3); in other words, 
IK/lloo < cll/Hoo, where || • ||oo denotes the uniform norm over f2. Since C n f and 
£* / are obtained by applying V to 2"*/, it follows immediately that 

|| £n || oo <c(logn) 2 and ||£* ||oo < c(logn) 2 , 

where n = (mod 3) and the uniform norm is taken over the triangle A. 

3.5.3. Fast Fourier transform. Comparing to the Hexagon-Hexagon case, the set 
up in the present subsection has at least one advantage if we consider the fast 
Fourier transform. The discrete Fourier transform of a function / periodic in H is 

W(t) = E where / k = (/, k )„ = -L £ /(^"^H 

For (•,•)„ in (|3.17p . we show that / k can be evaluated as in the classical discrete 
Fourier transform on a square. For this purpose, it is more convenient to use 
Cartesian coordinates. Let k = (^1,^2) corresponds to k. Then, by Theorem 12.21 

7k = (f,<M n = (f^ k )N = 4 E f(n- l Hj) f ? i ™- lk -1, 

ieA N 
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since <f> k (x) = e 2 *^ -1 *-* and B = nH- tr implies that cj) k (B- tr j) = e*™ -1 *^'. The 
homogeneous coordinates of Hj is (2j'i — 2j2 — ji, —ji —32), so that 

— 32 2 j 2 — jl — jl— J2 N-2t7TTt~ 1 fc-j 



A — __ \ ~* f /• 231-32 2 j 2 - ji -31-32 -v ' 
■ /k ~~ r,2 ■' ^ n ' n ' n / c 



3GA 



_ J_ \ ~" f( 2ji-32 232-31 -31-32 \ ' 
„2 ^ V n ' n ' n ^ 



232-31 —31—32 \„2i7TO 1 (fcl3'l+fc232) 



0<3li32<" 



This states that the discrete Fourier transform coincides, up to a reordering, with 
the classical discrete Fourier transform on a rectangle. Figure [9] shows the set A at 
and its reordering in rectangular coordinates. Similarly, recalling TL^ = {^k : k € 




Figure 9. The index set A^r (left) and its reordering (right). 



KjJ; the index set in rectangular coordinates can also be reordered, so that Hn 
becomes the product space in rectangular domain. In particular, this allows us to 
apply the classical FFT to evaluate /k- 

3.6. Other possibilities. There are other possible choices of lattices in our general 
frame of discrete Fourier analysis. For example, we can consider A = H~ tr and 
B = nH, for which the integral domain will be the hexagon in Figure 5. It 
is easy to see that the index sets A^r and Ajy in this case are and K n in the 
previous subsection, that is, their roles are interchanged. This case, however, does 
not seem to lead to interesting new result; the integral domain in the Stage 3 for 
the generalized cosine and sine functions will be the quadrilateral in Figure 7. 

One obvious question is if we can choose one lattice tiling R 2 with square or 
rhombus and choose the other lattice tiling E 2 with hexagon. The answer is neg- 
ative if we try to use regular hexagon, since the matrix H contains V3 and the 
requirement N = B tr A having all integer entries cannot be satisfied. We can, 
however, use other hexagon domains. For example, we can choose either 



Hi 



1 

-2 



Hi 



1 

-1 



Both lattices H1L 2 and B^i? tile R 2 . Their fundamental domains are depicted in 
Figure [101 The general result in Section 2.1 can be applied to develop a discrete 
Fourier analysis using either B\ or H.% and a lattice that tiles R 2 with either square 
or rhombus, since the requirement that N = B tr A has integer entries can be readily 
attained using, say A = H 1 or B 2 and B = I or R. Comparing to the regular 
hexagon, the hexagons in Figure 9 possess far less symmetry. The lack of symmetry 
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FIGURE 10. The fundamental domains of H\l? (left) and H 2 I? (right). 



means that we will not be able to carry the program outlined in Section 2.2 to Stage 
3 and Stage 4, whereas the results in Stage 1 and Stage 2 can be derived from the 
general theory straightforwardly. Hence, we will not pursuit the matter any further. 
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